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ABSTRACT 

A new generation of CMB experiments will soon make sensitive high reso- 
lution maps of the microwave sky. At angular scales less than ~10 arcminutes, 
most CMB anisotropies are generated at z < 1000, rather than at the surface of 
last scattering. Therefore, these maps potentially contain an enormous amount of 
information about the evolution of structure. Whereas spectral information can 
distinguish the thermal Sunyaev-Zeldovich (tSZ) effect from other anisotropies, 
the spectral form of anisotropies generated by the gravitational lensing and the 
kinetic Sunyaev-Zeldovich (kSZ) effects are identical. While spectrally identical, 
the statistical properties of these effects are different. We introduce a new real- 
space statistic, < 6(n) 3 6(rh) > c , and show that it is identically zero for weakly 
lensed primary anisotropies and, therefore, allows a direct measurement of the 
kSZ effect. Measuring this statistic can offer a new tool for studing the reion- 
ization epoch. Models with the same optical depth, but different reionization 
histories, can differ by more than a factor of 3 in the amplitude of the kSZ- 
generated non-Gaussian signal. 

Subject headings: cosmic microwave background — gravitational lensing — struc- 
ture formation 
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Introduction 



An upcoming generation of survey telescopes, such as ACT0, APExB and SPT§, will 
map the cosmic microwave background (CMB) with an unprecedented sensitivity and with 
resolution of ~1' . 

On degree scales, the main features of the CMB are fluctuations generated at the last 
scattering surface (z ~ 1100). However, on smaller angular scales, secondary effects pro- 
duced through the interaction between the CMB and the large-scale matter distribution of 
the universe dominate the microwave sky. Consequently, high resolution experiments will 
open a new window into different physical processes that affect the CMB photons as they 
travel from the recombination epoch to the present. In particular, these experiments are 
expected to provide invaluable information about the thermal Sunyaev-Zeldovich (tSZ), the 
kinetic Sunyaev-Zeldovich (kSZ) effects, and the weak gravitational lensing effects. 

The weak lensing effect appears to be a particularly promising probe of the growth of 
structure. Its imprints on the CMB sky have the potential of providing us with important 
information about the matter distribution over most of the cosmological history. So, besides 
testing our paradigm of structure formation, the measure ment of the lensing effec t can be 



very useful for constraining pa rameters like neutrino mass (IKaplinghat et a 



dark energy equation-of-state (IGoldberg fc Spergell Il999l ; IKaplinghat et al 



2003) and the 



2003J). Lensing 



measurem ents can also assist CM B polarization measurements in detecting primordial grav- 
ity waves (ISeljak fc Hiratal 12004 ) . 

There has been significant study of the effects of CMB lensing on the temperature and 
polarization power spectrum. Power spectrum analysis of t hese quantities can give us impor- 
tant information about the power spectrum of dark matter (jSeljakill996l ; IZaldarriaga Seljak 
19981 ). In addition, several works have been devoted to the development of methods for ex- 



tracting the projected matter field using t 
the CMB temperature and polarization ( 


le f 


ligher order correlations induced bv lensing on 


Hu 


2001; 


Hirata & Seliak 


2003a: 


Hu & Okamoto 


2002; 


Hirata & Seljak 


2003b: 


Okamoto & Hu 


2003 


)• 



Efforts to measure gravitational lensing may be limited by contamination produced by 
the tSZ and kSZ effects, and also by extragalactic poin t sources. Whereas these contaminants 
are thought to be unimportant for polarization maps (lHull2000[ ) , they will likely contaminate 
the temperature anisotropics at angular scales of ~ 10 arcminutes and below. The spectral 
dependence of the tSZ effect and the extragalactic point sources in principle provides the 



x http: / /www. hep. upenn.edu/act/act. html 
2 http: / /bolo. berkeley.edu/apexsz/ 
3 http://astro. uchicago.edu/spt/ 
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key for their removal from CMB temperature maps. However, the kSZ effect is spectrally in- 
distinguishable from the lensed primary CMB anisotropies, so a different method is required 
for removing this effect from CMB maps. 

A first method for reducing th e kSZ contamination of the lensing signal has already been 
proposed by lAmblard et al.l (12004 ). It uses the spatial correlation between the tSZ and kSZ 
effects to mask the main sources of kSZ anisotropies. Alternatively, in this work we suggest 
that the different statistical properties of the kinetic SZ and weak lensing effects, in particular 
the non-Gaussian features that they produce on the CMB anisotropies, can provide another 
tool for their separation. We introduce a new real-space statistic, < 9(n) 3 9(m) > c , that 
we show vanishes for weakly lensed primary CMB anisotropies. Consequently, this statistic 
would give us a clean measurement of the kSZ effect. This measurement is expected to be 
useful for quantifying and controlling the kSZ contamination of lensing field estimators, in 
addition to provide a new insight into the history of reionization and structure formation of 
the universe. 

This work is organized as follows: In chapter [2] we begin by reviewing the lensing ef- 
fect on the CMB temperature, and introduce the < 8(n) 3 8(rh) > c statistic, showing how 
it vanishes for weakly lensed primary anisotropies. In chapter [3] we review the kSZ effect 
and illustrate, using the Ostriker-Vishniac (OV) effect, how < 8 3 (n)6(m) > c can be used to 
determine the different physical parameters involved in our models for the kSZ effect. In 
chapter H] we present our conclusions. 



2. Weak Lensing of the CMB 

Weak lensing on the CMB temperature can be described as a mapping of the unlensed 
temperature anisotropies, 0(n), into the lensed ones, 9 L {n). (Note that we use 9{n) = 
T{n) /Tqmb — 1 for describing temperature anisotropies). The relation between these two 
maps is 

6 L (h) =9(n + 6h), (1) 

where the unitary vector n defines a direction on the sky and 5n represents the lensing de- 
flection angle due to the intervening mass distribution. 
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2.1. Weak lensing four-point function 

According to the standard inflationary scenario, the primordial density perturbations on 
the last scattering surface can be regarded as a nearly Gaussian random field. Consequently, 
the primary (unlensed) temperature anisotropies of the CMB will also have a Gaussian 
m-point probability distribution given by, 



P(6(n 1 ),...,6(n m ))d6(h 1 )...d8(h m ) = \ == ex P [ - \ ^X M ~ Vfa)l 

v /(27r) m det(M) lfr± x 

*dd{h 1 )...d9{h m ), (2) 



where m is an arbitrary positive integer and M, the covariance matrix, is such that My =< 
9(hi)6(hj) >. Therefore, the statistical nature of the unlensed anisotropies is completely 
specified by the two-point correlation function. More precisely, according to the Wick's 
theorem, any m-point correlation function for a Gaussian field will be given by, 

< t)( n , i . 6(n m ) >= <j ^P air s < 9 (^) e (f^) > ■ ■ ■< 9{n m ^)6{n m ) > for m even 

for m odd, 

where Epairs means that the sum must be performed over the (m — 1) ■ (m — 3) . . . 3 ■ 1 ways 
to make pairs. 

Therefore, the appearance of an additional term, usually called connected part, in the 
m-point function is a direct probe of non-Gaussianities that could have been induced by 
secondary effects like lensing or the kSZ effects. Since both effects are expected to have a 
null three-point function, in this work we will explore the potential of the connected part of 
the four-point function for accomplishing the separation of these two effects. 

In general, weak lensing produces a non-zero connected part for the four-point function 
of CMB temperature. It can be seen by Taylor expanding equation (JTJ), 

e L {n) w 6(h) + di0(h)(6n)i + ^0^9 (h) (5^(5^, (4) 

and calculating the first correction to the connected part of the four-point function due to 
lensing. Doing so one obtains, 

< 0L(ni)0L(7to)0x(n 3 )0L(n4) >c = \ E?,m )M =l < {^l)i{^m)j > X /gN 

< e(n p )di9(hi) >< 6(h q )dj6(hm) > &im P q, 



where the index "c" means "connected part", and Ai mpq is 1 if I, m, p, q are all different 
from each other and in all other cases. 
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2.2. The < 9(h) 3 9(m) > c statistic 

In this section, we introduce a new statistic, < 9(h) 3 9(m) > c , and show that it vanishes 
for the weakly lensed primary anisotropies. It is defined as 

< 9(h) 3 9(m) > c = < 9(h) 3 9(m) > -3 < 9(h) 2 >< 9(h)9(m) > . (6) 

Notice that, since it is a real-space statistic, it is not affected by any type of masking, and 
is straightforward to compute from temperature maps. 

As an illustration, the cancellation of this statistic for the weak lensing case can be 
easily seen from equation (jSj). If we set hi = h 2 = n 3 = n and ra 4 = m, we obtain, 

< 9 L (h) 3 9 L (m) > c = 6(< (SfyiiSh)] > + < (Sfyitfm^ >) < 6(m)di9(h) >< 9(n)dj9(h) > . (7) 

But the last factor of equation ((7j) is < 9{h)dj9{h) >= \ < dj9(n) 2 >= \dj < 9(h) 2 >= 0, 
provided that < 9(h) 2 > does not depend on h. This way we show that, to the lowest order 
in Sh, < 9L(h) 3 9i(hi) > c = 0. Eventhough equation ([3]) is a good approximation (since most 
of the lensing effect on the CMB is produced by linear density fluctuati ons, in which case 
\Sh\ is much smaller than the angular size of primary CMB anisotropies ( 1Seijaklll996l )). we 
provide below a more general demonstration. 

We know that the CMB anisotropies can be expanded into spherical harmonics, 

9 L (h) = 9(h + 6h) = J2 ai m Yr(n + 6h), (8) 

l in 

where J2im represents J2Zo Em=-r Then ' 

< 9 L (h) 3 9 L (m) > = <Y™ 1 (h + 5h)Y™ 2 (h + 5h)Y™ 3 (h + 5h)Y™ i (m + 5m) > 

iimi Z2?ri2 '3«^3 l4,m4 

X < dlimi^l2m2^'l3m3^l4m4 ^ • (9) 

Since we assume that the primary CMB anisotropies are a Gaussian random field, the last 
term in equation (jH]) can be decomposed into the sum of three terms corresponding to the 
three ways to make pairs. Thus, defining < ai m ay m i >= SfyS^^Ci where 5f v and 5^ m / 
correspond to Kronecker deltas and the bar over ai' m > represents the complex conjugate of 
ai' m i, it is easy to show that, 

< 9 L (h) 3 9 L (m) > = 3 GlGv < Y i m ^ + 5fl ) Y r{n + 5h)Y™' (h + Sn)Yf (m + 5m) > 

Ira I'm' 

oo oo I 

= CiCv < S Y ^ l[fi + 6fi ) Y i m ^ + 6fi ) 

1=0 l'=0 m=-l 

V 

x Y™' (n + ShftT 1 ' (m + 8m) > . (10) 
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But, making use of the addition theorem for spherical harmonics, 

1 - 2/ + 1 

V Y t m {n + 5h)Y^{m + 8m) = -^P z (cos( 7 )), (11) 

m=—l 

where Pi(x) corresponds to the Legendre polynomial and 7 is the angle between n + Sn and 
m + Sm, and also using that p(cos(0)) = 1, we can see that the sum Ylm=-i * n equation (fTUl) 
will be just a function of I and will not depend on the random field 5h, so < 6L(n) 3 8L{m) > 
can be separated into the product of two sums, 

00 21 4- 1 °° 2/' 4- 1 
< £ (n) 3 0i(m) >=3V^E C ''^ 1 - < fl'(cos( 7 )) > • (12) 

^— ' 47T 47T 

Z=0 J'=0 

But, using the definitions of aj m , Cj, and the addition theorem for spherical harmonics, it is 
easy to see that, 

00 9/ 4- 1 

< 8 L {n)8 L {m) >= ^ < P,(cos( 7 )) >, (13) 

so equation (fT2l) can be written, 

< e L (hf6 L (m) >= 3 < 6 L {nf >< e L (n)9 L (m) > . (14) 

This way we demonstrate that, 

< e L {nfe L {m) > c =< 9 L {nfe L {m) > -3 < 6 L (nf >< 6 L (n)6 L (m) >= 0. (15) 

It is important to keep in mind that the only approximation used in our calculations 
was to neglect any possible non-Gaussianity of the primary CMB anisotropies. We made no 
assumptions about the lensing deflection angle, Sn. Because of its cancellation, this statistic 
has the potential of giving us a direct measurement of the kSZ effect. 



3. The kinetic Sunyaev-Zeldovich effect 

The scattering of CMB photons by free electrons within galaxies and the intergalactic 
medium creates anisotropies on the CMB temperature. When these anisotropies are due 
to the thermal energy that the electrons transfer to the CMB photons, the effect is called 
thermal Sunyaev-Zeldovich (tSZ) effect. This effect modifies the thermal nature of the spec- 
trum of the CMB, reducing the number of low-energy photons and increasing the number of 
high-energy photons. At nearly 218 GHz the energy spectrum of the CMB remains almost 
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unchanged, giving rise to the so called thermal Sunyaev-Zeldovich "null", which is very im- 
portant for identifying this effect from other secondary CMB signals. 

On the other hand, if the free electrons have a bulk peculiar motion, additional temper- 
ature fluctuations are generated. This effect, known as the kinetic Sunyaev-Zeldovich (kSZ) 
effect, can be understood considering that in the rest frame of electrons the CMB is seen as 
hotter in one direction and colder in the opposite direction (Doppler effect). However, after 
their scattering with the electrons, part of the CMB photons are reemited isotropically in the 
rest frame of electrons generating a dipolar anisotropy in the rest frame of the CMB. This 
temperature aniso tropy, which does not depend on frequency, is given by (see, for example, 



Scannapiecol l2000l ) 



(«) = -\ 



7]rc 

dr}a(r})aTn e (nr], r])e~ T ^ ri ' > n ■ v(nr], 77), (16) 







where 77 represents the comoving distance along the path of CMB photons, which we will 
also use as a time coordinate, the subscript "rc" means "at recombination", a(rj) is the scale 
factor, v(r, rj) corresponds to the bulk velocity of electrons (note that throught this work 
bold letters will represent three-dimensional vectors, while two-dimensional vectors will be 
represented by symbols like v), c is the speed of light, n e (r, 77) is the electron density, a? is 
the cross section for Thomson scattering, and r(rj) = J" drj' aTn e (.v') a (.v') * s the optical depth 
to electron scattering at 77. Notice that in equation ffl6l) . like throughout this paper, we have 
assumed a flat universe. 

Since in this work we are interested in fluctuations at arcminutes scales, we can describe 
the kSZ effect by using the flat-sky approximation. Then, we will use x instead of n to refer 
to the angular position on the sky, where h will be assumed very close to the the z axis (i.e., 
h ~ (xi, X2, 1) where Xi <C 1). 

The kSZ effect is usually divided into two effects: the Ostriker-Vishniac (OV) effect 
and the non-linear kSZ effect. The OV effect is due to the presence of reionized matter 
within linearly evolving density fluctuations, while the non-linear kSZ effect is produced by 
electrons within non-linear structures. Given the highly predictive power of linear theory, 
the OV effect can give us accurate information about reionization, mainly in its initial phase. 
However, since in the low-redshift universe most of the reionized matter resides in high den- 
sity objects, a good understanding of both regimes is essential in order to have a complete 
picture of the kSZ effect. 
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3.1. The Ostriker-Vishniac effect 

In this work, we will focus only on the linear version of the kSZ effect, the OV effect. 
We will assume that the spatial distribution of free electrons is the same as the spatial 
distribution of dark matter, i. e., we will assume a homogeneously reionized universe and that 
dark matter traces baryons. Therefore, we will be interested in 5(r, if) = p(nr], r]) / p(r]) — 1 
where p represents the dark matter density. 

In the linear regime, the time evolution of the dark matter density perturbations is given 

by 

D(n) 

S(hr),r)) = 6(fir))j^, (17) 



where 5(nr]) = 5(nr], t] = 0), 



D{rj) = m v " / drj' 



2c J v ' a(r]')H(r]') 2 y ' 

is the growth function, and r/oo = r\[z = oo). 

It is useful to define the power spectrum P(k) of the density perturbations as 

< 6{k)5{k')* >= {27T) 3 P{k)5 D {k - k'), (19) 

where 5(k) and 5z?(k) correspond to the Fourier transform of 5(r) and the Dirac delta, 
respectively. 

In linear theory, the velocity field can be obtained from the continuity equation. Its 
Fourier transform is 

where D(i]) denotes the derivative of D(r]) with respect to time. 

From equation (|20|) we see that only Fourier modes with a non-zero line of sight com- 
ponent (i.e., k z in the flat-sky approximation) contribute to v z (r,i]). However, for a 
slow enough reionization process, the successive troughs and crests of these modes tend to 
cancel when projected along the z axis. Because of this, as we will see below, the OV effect 
arises from the modulation of the velocity field by the spatial fluctuations in the electron 
density. These fluctuations may be due to inhomogeneities in the baryon density and/or in 
the ionization fraction. Although the latter case, k nown as the patch y reionization effect, 



can give rise to important temperature anisotropies (ISantos et al.ll2003l ). we will concentrate 



in a homogeneously reionized universe where the OV effect is due only to fluctuations in the 
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baryon density. 

Then, in the flat-sky approximation, the temperature fluctuations can be written as 

d v{x) = — / dr)g(r))q z (xr],r),r)), (21) 

c Jo 

where g(rj) = a(r])aTn e (r])e~ T ^ is called the visibility function, and q(r, 77) = 6(r, r])v(r, rf). 
The Fourier transform of q(r, rf) is 

'^^Swf m ^-^ (22) 

From equation fl22|) . we can see that q z (k, rf) \k z =$ 7^ 0. This is very important because, given 
that the modes perpendicular to z do not cancel out when projected along the line of sight, 
in a slowly varying reionization process they dominate the OV effect. 
Finally, it can be shown that the Fourier transform of 9{x) is 

MO = - r d^M [ ^e-^q z (-, K rj), (23) 
J a] 2 J 2tt 7? 

where l/rj and k z are the components perpendicular and parallel to the z axis of k. 



3.2. The < 0(xi) 3 6(x 2 ) > c statistic for the OV effect 



For the calculation of the < 6(xi) 3 6(x*2) > c statistic for the Ostriker-Vishniac effect, 
< ov (xi) 3 9 ov (x~2) > c , we must compute the trispectrum T m {li, l 2 , 13, which is defined by 



< Qov(,h)9ov(,h)9ov(h)9ov(h) >c 



(2n) 2 T ov {h,l 2 ,l 3 ,k)5(h + l 2 + l 3 + U) 



(24) 



where 9 ov (l) is given by equation ((231) . In order to do so, we will use the Limber's ap- 
proximation ( Kaiserl Il992l ; ICastrol 12004 ). It consists in considering 9 ov (x) as the sum of 
the contributions A8 ov (x) from many shells centered at rj and of width Arj which is much 
larger than the relevant wavelengths. These conditions can be expressed as Ai]/i] rc 1 and 
Arj/r] ^> ip where ip represents the angular scale of interest. Then, according to equation 
the Fourier transform of A6 nv (x) is 



A9 ov {l) 



'/o- 



A?7 



'/() - 



A>7 



cr}'- 



2tt 



V 



(25) 



If g(rj) and g 2 (-, k z , rj) vary slowly with the comoving distance r], they can be considered 



constant within each shell. Then, one can obtain that 



A9 0V {1) 



g{r)o)Ar] f dk 



evil 



2tt 



^k z , Vo ) Jo (k z ^ 



2 J 



(26) 
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where jo(x) = sin(x)/x. 

In the Limber's approximation, the fields defined within different shells are, by assump- 
tion, statistically independent. Therefore, the trispectrum of the OV effect is the sum of the 
independent trispectra generated in all the shells. 

For the calculation of these trispectra, we need to compute < §*(—, k z i, r))...q z (^, k z ±, rj) >. 
In order to do so, we use the fact that, in the linear regime, 8(r) can be considered as a 
Gaussian random field, and use the Wick's theorem. 

The presence of jo{k z ^-) in equation (l2"B"j) implies that the contribution to the ef- 
fect comes mainly from k z < 1/Arj. In addition, the condition At]/i] ^> (p means that 
1/t}^> I/A77. Therefore, as stated above, the main contribution to the OV effect comes from 
Fourier modes of q z (k, if) nearly perpendicular to the line of sight. Then, we can approximate 
q z (-, k z , rj) ~ q z ( -, 0, rf) and perform the integration over k z in equation f[2"6"j) in a very easy 
way. 

Finally, approximating the sum of the trispectra from all the shells to an integral, one 
can obtain that 



T ov (h, h, h, h) — q 

o 

where 

f(h,k 2 ,k 3 ,h) = 



no 



dr\- 



D 2 



Tj Tj Tj Tj 



(27) 



dk 



(27T) 



P(a)P(b)P(c)P(d)k 4 z 



i,j,k=l 



X 



(28) 

a = |k|, b = |k — fcfel, c = [k + fcj and d = |k + fej + kj\. Notice that equation (1271) differs 
from equation (60) of ICastrol (120041 ) by the combinatorial factor 1/4!. 
Now, if we define 



< M^i) 3 M^2) > c = 



dl 



ft 



/■^-^c e l\i), 



we have 



c e : e (i) 



ft 



di{ 



d% 



(2vr) 2 J (2tt) 



■T ov (h, h — hi ~h — 1,1)- 



(29) 



(30) 



So, considering equation (l2"7j) . we obtain 



d g{r]Y faD(rj)D(rj)Y f dh dk 2 dk 



D 2 



(2vr) 2 (2tt) 2 (2tt) 3 



(31) 
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where a! = |k|, b' = |k — l/r]\, d = |k — kx\, and d! — |k — k<i 



3.3. Discussion 

The amplitude and shape of C^ e (l) depends on the reionization history of the universe 
(contained in the visibility function g{f])). There is not yet certainty about the actual shape 
of g(i]); in fact, the aim of this work is to contribute to answering this question. However, 
in order to show the kind of information that can be obtained from the < 8(n) 3 8(m) > c 
statistic, in this section we calculate C^f{l) making use of a very simple model for g(rj). 

In order to do so, we use the flat-ACDM cosmo logical model with the cosmo logical 
parameters given by (ft m , fib, h, n. s , r, er a) = (0-24, 0-04, 0.7 3, 0.95, 0.09, 0.74), in agreement 



with the three years WMAP results (ISpergel et aJj l2006l) . For the transfer function we 



make use of the fitting form given by iBardeen et al.l (I1986I ). We assume a mass fraction of 



hydrogen and helium of 75% and 25%, respectively, with simultaneous hydrogen and helium 
reionizations. 

For the reionization history, we assume a very simple two-stage model. In this model 
the universe is partially reionized at redshift z r i, with an ionization fraction of x° e . After 
that, at z < 7, the universe is completely ionized, i.e., 

Z > Z r i 

x e (z) = { x° e z ri > z>7 (32) 

1 z < 7. 

In figure ([1]), we show the calculation of Cf v e (l) for four different reionization scenar- 
ios. The solid, long-dashed, dashed, and dotted lines, represent reionization histories with 
(x° e ,z ri ) = (0.8, 13.3), (0.6, 15.1), (0.4, 18.4), and (0.2,27), respectively. Notice that all these 
scenarios have the same optical depth to Thompson scattering at the last scattering surface, 
r = 0.09 and the same cr 8 = 0.74. 

We can see that, at angular scales relevant for the upcoming high resolution CMB 
experiments (10 3 < / < 10 4 ), C^ e (l) varies by approximately a factor of 3 between the 
comparatively early [z r i = 27) and late (z« = 13.3) reionization scenarios. This result sug- 
gests that, eventhough < 9(n) 3 9(m) > c also depends on other cosmological parameters like 
t and as, it would be sensitive enough to provide us with new information about the history 
of reionization of the universe when complemented with other independent measurements. 
However, before arriving to any conclusion about the relationship between C^ e (l) and the 
physical parameters involved in reionization, our calculation must include more realistic 
models for the reionization process, and also consider the non-linear behaviour of the kSZ 
effect. 
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Fig. 1. — Z 4 C^ 9 (/) for four two-stage reionization scenarios. The solid, long-dashed, dashed, 
and dotted lines, represent reionization histories with (x®,z r i) = (0.8,13.3), (0.6,15.1), 
(0.4,18.4), and (0.2,27), respectively. In the four cases we have assumed r = 0.09 and 
a 8 = 0.74. 
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4. Conclusions 

We propose non-Gaussianity to separate the kSZ effect from the lensing effect in CMB 
maps. Specifically, we show that a particular shape of the trispectrum, < 9 3 (n)9(rh) > c , 
cancels for weakly lensed primary CMB anisotropies providing a very good tool for measur- 
ing the kSZ effect. 

Measurements of the amplitude of the kSZ signal complement measurements of the op- 
tical depth and provide new insights into the reionization history of the universe. We show 
that reionization histories that produce similar low I < EE > signal have very different kSZ 
signatures. In addition to yield a new insight into the history of reionization and structure 
formation of the universe, < 9 3 (h)9(m) > c could become very useful for quantifying and 
controlling the kSZ contamination of lensing statistics. 

While our calculation is based on a linear theory estimate and needs further investiga- 
tion with full numerical simulations of reionization, we suspect that they will confirm our 
basic conclusion that kSZ measurements will shed more light on the dark ages. 
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